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ABSTRACT 


Atmospheric  turbulence  degrades  a  coherent  laser  beam  uihen  it  propagates  through  the 
atmosphere.  Measurements  of  the  distribution  of  atmospheric  turbulence  provide  insight  into  the 
underlying  mechanisms  that  produce  optical  turbulence  and  suggest  possible  means  to  overcome  or 
circumvent  the  effects  of  such  turbulence.  A  variety  of  acoustic,  optical  and  thermal  probe 
instruments  provide  measurements  of  atmospheric  turbulence.  Of  these,  the  acoustic  echosounder  can 
measure  atmospheric  density  and  velocity  irregularities.  During  the  course  of  previous  work, 
questions  arose  concerning  the  calibration  of  the  NPS  echosounder.  The  echosounder  appears  to 
detect  a  higher  level  of  atmospheric  turbulence  than  do  other  instruments  used  at  the  same  altitude. 
This  resulting  overestimation  of  atmospheric  turbulence  could  significantly  influence  programs  such 
as  the  Advanced  Electro-Optic  Site  (AEOS),  a  proposed  4  m  telescope  to  be  built  in  Hawaii. 

This  thesis  will  attempt  to  identify  key  components  contributing  to  the  sounder  calibration 
including  the  transducer  transmitting  and  receiving  efficiencies  as  well  as  their  dqiendence  on 
pressure  and  atmospheric  density.  Additionally  a  modified  acoustic  echosounder  equation  will  be 
developed  which  more  accurately  profiles  the  atmospheric  turbulence  measured  by  such  a  sounder. 
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I.  INTRODUCTION 


Turbulence  degrades  a  coherent  laser  beam  when  it  propagates  through 
the  atmosphere.  Atmospheric  irregularities  randomize  the  amplitudes  and 
phases  of  the  electromagnetic  waves  as  they  propagate  through  these  turbulent 
regions.  The  result  is  that  an  initially  focused  beam  is  defocused,  which 
reduces  the  irradiance  delivered  to  the  target. 

Measurements  of  the  distribution  of  atmospheric  turbulence  provide 
insight  into  the  underlying  mechanisms  that  produce  optical  turbulence  and 
suggest  possible  means  to  overcome  or  circumvent  the  effects  of  such 
turbulence. 

A  variety  of  acoustic,  optical  and  thermal  probe  instruments  provide 
measurements  of  atmospheric  turbulence.  Of  these,  the  acoustic  echosounder 
can  measure  atmospheric  density  and  velocity  irregularities  resulting  from  air 
currents,  temperature  inversions,  humidity  variations,  mechanical  turbulence 
and  other  causes. 

Walters  has  probed  the  atmosphere,  quantitatively  characterizing  its 
turbulence,  over  the  last  decade.  He  has  made  extensive  measurements  in  the 
mountainous  area  in  the  vicinity  of  White  Sands,  New  Mexico  [Ref.  1]. 

Working  with  Walters,  Weingartner  [Ref.  2]  and  Wroblewski  [Ref.  3]  built 
and  tested  a  25  element  5  kHz  planar  array  acoustic  echosounder.  Moxcey 
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[Ref.  4]  later  modified  this  design  by  placing  the  drivers  closer  together  into  a 
19  element,  hexagonal,  close-packed  array.  This  modification  increased  the 
acoustic  power  in  the  main  lobe  and  suppressed  the  side  lobes.  The  result  was 
an  extremely  effective  low  altitude  echosounder  that  was  less  expensive,  more 
portable  and  had  better  side  lobe  reduction  than  previous  monostatic  sounders. 

During  the  course  of  the  previous  work,  questions  arose  concerning  the 
calibration  of  the  echosounder.  The  echosounder  appears  to  detect  a  higher 
level  (by  a  factor  of  almost  four)  of  atmospheric  turbulence  than  do  other 
instruments  (e.g.,  high  speed  temperature  probes)  used  at  the  same  altitude. 

The  resulting  overestimation  of  atmospheric  turbulence  could  significantly 
influence  programs  such  as  the  Advanced  Electro-Optic  Site  (AEOS),  a 
proposed  4  m  telescope  to  be  built  on  Mount  Haleakala,  Maui,  Hawaii.  For 
AEOS  it  is  crucial  to  find  the  best  design  compromise,  to  reduce  atmospheric 
turbulence,  either  naturally  or  by  building  a  telescope  support  pedestal  above 
the  turbulent  layer.  An  overestimation  of  turbulence  would  require  a  pedestal 
higher  than  is  actually  necessary  to  get  above  the  turbulent  layer,  consequently 
wasting  money. 

This  thesis  will  attempt  to  Identify  key  components  contributing  to  the 
acoustic  sounder  calibration,  including  the  transducer  transmitting  and  receiving 
efficiencies,  as  well  as  their  dependence  on  pressure  and  atmospheric  density. 
Measurements  made  in  the  NFS  anechoic  chamber  will  be  compared  with  field 
calibration  measurements  made  at  the  Starfire  Optical  Range  against  three 


2 


differential-temperature  probes  mounted  on  a  tower.  A  modified  acoustic 
echosounder  equation  will  be  developed  that  utilizes  transducer  calibration 
parameters  measured  in  the  laboratory. 
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il.  BACKGROUND 


A.  ATMOSPHERIC  TURBULENCE 

Random  phase  perturbations  occur  as  optical  waves  propagate  through 
the  atmosphere.  These  small  scale  phase  variations  cause  the  lirht  beam  to 
diverge,  reducing  the  irradiance  on  a  target.  Atmospheric  density  fluctuations 
carried  by  the  turbulent  velocity  field  cause  these  optical  phase  perturbations 
by  affecting  the  local  index  of  refraction.  Temperature  fluctuations  create  such 
density  variations,  which  can  occur  over  scales  anywhere  from  tens  of  meters 
in  size  down  to  sub-centimeter  ranges. 

These  phenomena  degrade  the  resolution  of  both  earth-bound  observation 
facilities  and  laser  platforms.  The  following  two  parameters  characterize  this 
degradation:  the  refractive  turbulence  structure  parameter,  and  the  spatial 
coherence  length  of  the  atmosphere,  r„.  These  parameters  are  most  variable 
and  inhomogeneous  at  the  interface  between  stratified  atmospheric  layers.  This 
is  particularly  relevant  near  the  earth's  surface. 

1 .  Refractive  Turbulence  Structure  Parameter 

The  refractive  turbulence  structure  parameter,  C„^,  is  the  mean- 
squared  statistical  average  of  the  difference  between  the  indices  of  refraction 
at  two  points,  scaled  by  the  points'  separation.  Tatarski  [Ref.  5}  includes  a 
detailed  discussion  of  this  parameter,  defined  as 
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(1) 


^2_  <ini-Ih)^> 

«  ■  r  2/3 
-^12 

where  <  >  signifies  an  ensemble  average,  n^  and  nj  are  the  indices  of 
refraction  at  points  1  and  2  respectively,  and  r^j 's  the  separation  between 
points  1  and  2  in  the  atmosphere.  is  a  measure  of  the  local  variability  in 
refractive  index,  and  it  provides  the  starting  point  in  analyzing  the  resulting 
optical  turbulence.  It  is  important  to  note  that  the  presence  of  significant 
velocity  turbulence  does  not  imply  an  appreciable  degree  of  optical  turbulence, 
and  vice  versa  [Ref.  6].  Unfortunately,  is  a  difficult  parameter  to  measure 
directly.  It  is  easier  to  make  direct  measurements  of  another  atmospheric 
structure  parameter,  Ct^. 


2.  Temperature  Structure  Parameter 

The  temperature  structure  parameter,  C/,  is  the  analogous  measure 
of  temperature  fluctuations  between  two  points  in  space 


C  2  = 
'-T 


r  2/3 
■‘■12 


where  T^  and  Tj  are  the  temperatures  at  points  1  and  2  respectively  and  r^j  is 


defined  as  in  Equation  (1).  As  will  be  discussed  in  the  following  section,  an 


acoustic  echosounder  measures  this  parameter  indirectly.  For  optical 
wavelengths,  Tataraki  [Ref.  5]  developed  the  following  relationship  between 
and  Ct^: 
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(3) 


p  2  _  /  "79x10  ^  P  \  2  2 

where  P  is  the  atmospheric  pressure  in  millibars,  and  T  is  the  temperature  in 
kelvin.  Equation  (3)  ignores  a  contributing  factor  that  depends  on  the  water 
vapor  concentration,  that  is  important  quantifying  the  effects  of  turbulence  on 
the  propagation  of  microwave  radiation. 

3.  Spatial  Coherence  Length 

As  previously  mentioned,  the  spatial  coherence  length,  r„,  introduced 
by  Fried  [Ref.  7],  is  an  important  optical  turbulence  parameter.  It  is  a  measure 
of  the  spatial  correlation  of  the  electric  fields  of  the  optical  beam,  measured  in 
a  plane  perpendicular  to  the  direction  of  propagation,  to  is  a  function  of  the 
magnitude  of  the  integrated  optical  turbulence  and  depends  on  the  integral  of 
over  the  propagation  path  in  question.  For  a  plane  wave, 

=  2.1  [1.46 /c2  cfz] -3/5  ,  (4) 

where  k  is  the  wave  number  (2/rM)  and  L  is  the  optical  path  length.  Coherence 
lengths  for  vertical  beam  propagation  vary  widely  between  observation  sites. 
Coherence  length  values  can  range  from  a  few  centimeters  at  a  poor  site  to 
around  30  cm  at  a  world-class  site. 
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B.  THE  ACOUSTIC  ECHOSOUNDER 

Acoustic  echosounders  transmit  a  pulse  of  acoustic  energy  to  detect 
atmospheric  density  variations  in  much  the  same  way  as  a  shipborne  sonar 
detects  the  presence  of  a  submarine.  The  density  structure  of  the  atmosphere 
scatters  a  portion  of  the  transmitted  energy  back  to  the  echosounder.  The 
power  returned  from  these  turbulent  structures,  Ph,  is  proportional  to  the  power 
transmitted,  P^.  The  echosonde  equation,  also  known  as  the  radar  equation, 
has  been  adapted  by  Neff  [Ref.  8],  based  upon  work  by  Tatarski  [Ref.  5],  Little 
[Ref.  91  and  Hall  [Ref.  10],  to  the  operation  of  an  echosounder: 


-^=[p^rj  [e-2«*l  [-^l  [-^G] 


T-“T- 


(5) 


with. 


•  Pr/Er  is  the  received  acoustic  power,  where  Pr  is  the  measured  electrical 
power  and  Er  is  the  efficiency  of  conversion  from  received  acoustic  power 
to  electrical  power. 

•  P-tEt  is  the  radiated  acoustic  power,  where  P^  is  the  electrical  power 
applied  to  the  transducer  and  E^  is  the  efficiency  of  conversion  of 
electrical  driving  power  to  radiated  acoustic  power. 

•  e'^'^  is  the  round  trip  power  loss  resulting  from  attenuation  by  air,  where 
a  is  the  average  attenuation  (m  ')  to  the  scattering  volume  at  range  R  (m). 

•  acoustic  scattering  cross-section  per  unit  volume;  that  is,  the 
fraction  of  incident  power  backscattered  per  unit  distance  into  a  unit  solid 
angle  at  frequency  f  (Hz). 

•  cr/2  is  the  maximum  effective  scattering  volume  thickness,  where  c  is  the 
local  speed  of  sound  (ms  and  r  is  the  pulse  length  (s). 
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•  AG/R^  is  the  solid  angle  subtended  by  the  antenna  aperture  A  (m^)  from 
the  scattering  volume,  modified  by  an  effective  aperture  factor  G,  arising 
from  the  antenna's  directivity  and  geometry. 

As  mentioned  previously,  the  acoustic  echosounder  measures  atmospheric 
turbulence,  values  via  the  scattering  cross-section,  a.  Tatarski  [Ref.  5] 
shows  that  the  backscattering  cross-section  is  proportional  to  through  the 
following  relationship, 

a{R,f)  =0.0039ic^/3  (£r  )  (5) 

where, 

•  T  is  the  absolute  temperature  (K), 

•  k  Is  the  acoustic  wave  number  {m"’), 

•  vt  is  the  acoustic  wavelength  (m). 

This  is  the  case  only  for  backscattered  acoustic  energy,  that  is,  energy  which 
is  scattered  through  an  angle  of  180  degrees.  For  other  angles,  the  amount  of 
energy  scattered  depends  upon  the  velocity  structure  parameter  of  the 
atmosphere  [Ref.  51.  Combining  Equations  (5)  and  (6),  and  solving  for  Ct^ 
gives. 


c  2  =  [. 


^  0.0039  Ej.EgAG^ 


r  Pn 

lci/3  ‘  CT  Pj 


(7) 


Several  of  the  factors  appearing  in  Equation  (7),  such  as  P^,  E^,  E^  and  G  are 
difficult  to  determine  accurately.  Therefore,  a  certain  degree  of  error  exists 
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with  an  acoustic  sounder  that  must  be  dealt  with  when  calibrating  such  a 
device.  One  calibration  technique  involves  operation  of  an  acoustic 
echosounder  in  the  vicinity  of  a  tower  that  has  high  speed  temperature  probes 
mounted  at  known  heights.  The  measurements  from  the  temperature  probes 
are  compared  to  the  echosounder  data.  Although,  straightforward,  this 
procedure  does  not  allow  for  extrapolation  of  the  calibration  to  other  locations, 
particularly  when  significant  atmospheric  density  changes  occur.  Reference  1 1 
discusses  echosounder  calibration  in  detail. 

The  heart  of  the  NFS  echosounder  Is  an  array  of  1 9  Motorola  KSN  1 005A 
piezoelectric  speakers  mounted  in  a  hexagonal  shaped  array.  These  speakers 
were  chosen  for  their  resonance  properties  and  high  transmission  and  receiving 
efficiencies.  The  complete  echosounder  system  is  shown  in  Figure  1 . 
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III.  DISCUSSION 

A.  PREVIOUS  WORK 

1 .  Transducer  Efficiency 

The  acoustic  echosounder  is  an  excellent  tool  for  probing  the  lower 
troposphere,  however  it  has  a  limited  range  ( ^  200  m)  that  depends  on  the 
acoustic  frequency  used  and  on  atmospheric  attenuation.  The  system 
estimates  atmospheric  turbulence  by  measuring  the  absolute  acoustic  energy 
scattered  by  the  atmosphere.  The  transducer  transmitting  and  receiving 
efficiencies  must  be  known  in  order  to  measure  the  returned  power. 
Transducer  efficiency  can  typically  range  anywhere  from  a  few  percent  [Ref. 
8]  to  25  percent  [Ref.  121. 

The  echosounder  operation  is  frequency-dependent.  Its  ability  to 
measure  the  spatial  resolution  of  temperature  structures  increases  at  higher 
operating  frequencies,  but  atmospheric  attenuation  increases  exponentially  with 
frequency  as  well.  The  Motorola  KSN  1005 A  speakers  used  in  the  hexagonal 
array  have  a  maximum  response  at  a  resonant  frequency  of  5000  Hz.  Although 
a  lower  frequency  would  reduce  atmospheric  attenuation,  5000  Hz  optimizes 
the  speaker  efficiency  and  spatial  resolution  [Ref  2]. 

Weingartner,  under  Waiters'  supervision,  assembled  the  first 
operational  acoustic  echosounder  system.  He  measured  the  transducer 
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receiving  efficiency,  Er,  and  found  it  to  be  50  percent  [Ref.  2].  This  was  much 
higher  than  the  efficiency  predicted  by  References  8  and  1 2.  He  was  not  able 
to  measure  the  transmitting  efficiency,  E^,  so  he  assumed  that  the  values  for 
Er  and  Et  were  equal,  with  value  0.5. 

2.  Effective  Aperture  Factor 

Weingartner's  system  employed  a  5  x  5  (25  element)  square  array 
of  transducers.  He  used  calculations  that  Hall  and  Wescott  [Ref.  10]  made  for 
the  beam-shape  compensation  factor.  Their  work  was  based  on  the  theoretical 
development  by  Probert-Jones  [Ref.  13].  This  value,  0.40,  is  the  same  as  the 
effective-aperture  factor,  G,  in  Equation  (7). 

Moxcey  [Ref.  4]  modified  the  array  by  reducing  the  number  of 
speakers  from  25  to  19  and  packing  these  transducers  closer  together  into  a 
hexagonal  configuration.  This  reduced  the  side  lobes  and  increased  the  power 
in  the  main  lobe.  Approximating  the  aperture  area  of  the  array  to  be  equivalent 
to  19  times  the  aperture  area  of  a  single  speaker  having  a  diameter  of  7.620 
centimeters,  one  gets  a  value  of  0.0866  m^  for  the  area  of  the  array.  This  is 
the  antenna  area.  A,  in  Equation  (7). 

Equation  (7)  uses  the  above  values  for  Er,  Et,  G  and  A.  Combining 
these  values  with  the  constants  and  simplifying  we  get, 

C/ =  59217  [-^]  [-^ly]  ,  (8) 

for  the  relationship  between  and  the  echosounder  system. 
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B.  FIELD  MEASUREMENTS 


1 .  Advanced  Electro-Optic  Site  (AEOS) 

The  Air  Force  is  planning  to  build  a  4  m  telescope  on  the  top  of 
Mount  Haleakala  on  the  island  of  Maui  in  Hawaii.  Turbulence  will  interfere  with 
the  "seeing  ability"  of  the  telescope.  Previous  data  indicate  that  a  pedestal  can 
raise  the  telescope  above  the  majority  of  the  low  level  turbulence.  This  appears 
to  be  feasible  since  the  majority  of  the  turbulence  is  confined  to  a  ~  20  m  level 
layer  which  hugs  the  contours  of  the  mountain.  Two  questions  are:  how  thick 
is  this  ground  layer  and  how  much  turbulence  does  it  contain?  Also,  how  high 
does  the  pedestal  have  to  be  to  get  above  this  ground  layer? 

Measurements  with  an  acoustic  echosounder  were  made  on  three 
separate  occasions  to  determine  the  extent  of  the  turbulence  in  the  ground 
layer.  The  first  measurements  were  made  by  the  Environmental  Studies 
Section  of  the  Rome  Air  Development  Center  jRADC)  in  August  1 974  (Ref.  1 4]. 
Their  echosounder  did  not  find  any  significant  turbulence,  of  sufficient  strength 
to  degrade  seeing,  within  the  range  100-1000  ft.  Their  data  for  the  range  0- 
100  ft  was  inconclusive. 

The  second  set  of  data  reported  by  Mattingly  [Ref.  1 1  ]  were  recorded 
during  the  period  26  June  -  2  July  1991 .  He  found  an  intensely  turbulent  layer 
within  the  first  80  meters  above  AEOS.  He  calculated  a  30-40%  improvement 
in  both  daytime  and  nighttime  seeing  quality  for  a  1 5-25  m  tower. 


The  third  set  of  data  reported  by  Gast  [Ref.  1 51  were  recorded  during 
the  period  3-12  March  1992.  His  measurements  indicated  significant 
turbulence  present  in  the  region  from  10-100  m  above  AEOS,  although  this 
turbulence  was  not  as  intense  as  that  found  by  Mattingly  [Ref.  11].  Gast  found 
that  the  echosounder  was  reporting  turbulence  at  a  level  approximately  four 
times  larger  than  would  be  expected.  This  was  based  upon  the  maximum  value 
that  Walters  measured  for  the  log  of  C„^  (-13.0),  in  the  desert,  at  noon  when 
atmospheric  turbulence  is  greatest  [Ref  1  ].  Figure  2  shows  that  the  AEOS  data 
exceeded  this  maximum  value  for  over  ten  hours  20  m  above  the  ground. 
However,  problems  in  the  microthermal  probe  electronics  used  for  calibration 
make  these  comparisons  suspect  and  inconclusive. 
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Local  Time  (Hours) 


Figure  2.  data  from  AEOS  which  shows  that  unmodified  measurements 
exceed  maximum  expected  values  (-13.0)  for  a  desert  location  at  high  noon, 
the  most  turbulent  time  of  day. 
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2.  Starfire  Optical  Range  (SOR) 


Measurements  were  conducted  at  the  Starfire  Optical  Range,  Kirkland 
Air  Force  Base,  Albuquerque,  New  Mexico  on  two  occasions.  The  first  set  of 
measurements  took  place  from  28  through  30  April  1992.  The  calibration  for 
these  measurements,  at  11 .5  and  40  m  altitude,  was  found  to  be  invalid  when 
it  was  discovered  that  an  RMS  chip  inside  the  reference  probes  had  a  floating 
connection  that  introduced  a  substrate  bias  error.  The  reference  probe  used  at 
20.5  m  altitude  operated  correctly,  so  those  data  are  not  suspect.  The  floating 
ground  problem  was  discovered  and  corrected  in  July  1992. 

The  second  set  of  measurements  occurred  during  the  period  5-9 
September  1992.  Reference  probes  were  positioned  at  1 1.5,  20.5  and  40  m 
altitudes.  Figures  3  through  5  show  that  all  probes  operated  correctly.  This 
was  the  first  time  that  the  acoustic  echosounder  was  compared  directly  with 
a  credible  independent  calibration.  However,  a  factor  of  four  discrepancy 
between  the  two  data  sets  was  evident.  Figures  3  through  5  show  that  the 
echosounder  values  closely  resemble  the  tower  mounted  high  speed 
temperature  probes,  after  the  echosounder  data  has  been  divided  by  four. 
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5-6  Sep  92  Calibration  Run 


Local  Time  (Hours) 


- Tower  1  1.5  m  - Sounder  12.5  m 


Figure  3.  Echosounder  data,  divided  by  four,  resembles  data  collected 
simultaneously  by  tower-mounted,  high  speed  temperature  probes  at  a 
similar  height. 
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5-6  Sep  92  Calibration  Run 


Locol  Time  (Hours) 


- Tower  20.5  m  - Sounder  20  m 


Figure  4.  Echosounder  data,  divided  by  four,  closely  resembles  data 
collected  simultaneously  by  tower-mounted,  high  speed  temperature  probes 
at  a  similar  height. 
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5-6  Sep  92  Calibration  Run 


Local  Time  (Hours) 


- Tower  40  m  - Sounder  40  m 


Figure  5.  Echosounder  data,  divided  by  four,  closely  resembles  data 
collected  by  tower-mounted,  high  speed  temperature  probes  at  a  same 
height. 
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C.  ECHOSOUNDER  OPERATION 

The  echosounder  can  only  measure  atmospheric  turbulence  accurately 
when  the  system  performance  factors  (i.e.,  Er,  Ej,  Pj,  G)  are  known  and 
incorporated  into  Equation  (7).  These  factors  will  be  measured/calculated  to 
determine  if  they  are  responsible  for  the  echosounder  reporting  turbulence  at 
a  level  four  times  larger  than  would  be  expected.  Additionally,  the  transducer 
efficiency  dependence  on  pressure  was  investigated  in  order  to  use  the  SOR 
echosounder  calibration  to  evaluate  data  at  a  different  pressure  e.g.,  AEOS. 

1.  Receiving  Efficiency 

The  receiving  efficiency  was  measured  in  the  NPS  anechoic  chamber, 
using  a  HP3314A  Function  Generator,  a  HP  3561 A  Dynamic  Signal  Analyzer, 
a  one  inch  reference  speaker  and  the  19  element  hexagonal  array.  Figure  6 
shows  the  experimental  equipment  arrangement.  The  function  generator  with 
a  50  ohm  output  impedance  provided  an  open  circuit,  1.0  volt,  5.0  kHz 
sinusoidal  signal  to  a  one  inch  dome  driver.  The  array  was  aligned  along  the 
centroid  of  the  main  lobe  of  the  reference  speaker  signal.  The  array  received 
the  transmitted  acoustic  signal  and  the  HP  3561 A  Signal  Analyzer  measured 
the  RMS  amplitude  of  the  signal. 

The  array  was  then  replaced  with  a  1/4  inch.  Cartridge  Type  4135, 
Bruel  &  Kjaer  (B&K),  calibrated  condenser  microphone,  with  Type  2633 
preamplifier.  The  microphone  was  calibrated  using  a  B&K  pistonphone,  prior 
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1 ■  Driver 

r  n  II  t  n  u  t  ^  v«A - 1 


Figure  6.  Receiving  efficiency  measurement  equipment  set-up  inside  the 
anechoic  chamber. 
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to  use.  The  pistonphone  has  a  chamber  of  accurately  known  volume  in  which 
a  piston  of  known  area  oscillates  with  known  displacement  amplitude  [Ref  1 6]. 
The  microphone  calibration  coefficient,  Cpj«o„,  was  the  RMS  signal  voltage 
divided  by  the  RMS  pressure  generated  by  the  pistonphone  calibration  source. 
The  measurements  made  with  the  microphone  were  then  divided  by  this 
coefficient  to  scale  the  data. 

The  microphone  determined  the  total  power  available  for  reception 
by  the  acoustic  array,  i.e.,  the  total  power  actually  ensonified  within  the  solid 
angle  of  the  acoustic  array.  An  HP  3561 A  digital  spectrum  analyzer  made  the 
microphone,  RMS  voltage  measurements.  With  the  small  one  inch  dome 
speaker  providing  a  reference  acoustic  signal,  the  calibrated  microphone 
measured  the  sound  pressure  field  with  the  echosounder  array  removed.  After 
placing  the  array  in  the  sound  field,  the  receive  sensitivity  was  the  open  circuit 
voltage  divided  by  the  incident  pressure. 

2.  Transmitting  Efficiency 

The  transmitting  efficiency  was  harder  to  measure.  Again,  this 
measurement  was  made  inside  the  anechoic  chamber.  The  HP3314A  Function 
Generator  generated  a  1 .0  Vp^,  5.0  kHz,  sinusoidal  signal  with  a  50  Q  output 
impedance  which  was  connected  to  the  hexagonal  array  i.e.,  the  array  was 
used  to  transmit  the  signal  this  time.  Figure  7  shows  the  experimental 
equipment  arrangement.  A  1/4  inch.  Cartridge  Type  4135,  B&K  calibrated 
condenser  microphone,  positioned  on  a  thin  rod,  received  the  signal  and  passed 
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Figure  7.  Transmitting  efficiency  measurement  equipment  set-up  inside  the 
anechoic  chamber. 


it  to  the  signal  analyzer.  A  signal  amplitude  measurement  was  made  at  the 
centroid  of  the  main  lobe  and  then  again  at  5  cm  increments  out  to  the  first 
node. 

These  voltage  measurements  (VJ  were  first  converted  to  pressure 
using  the  conversion  coefficient  {C^„)  from  the  pistonphone 

V 

p  =  ^  rioi 

^presaure  ^ 

^piston 

where  is  the  RMS  pressure  in  pascal.  This  pressure  was  then  converted 

to  a  power  density  (I)  using  the  relationship 


^praBBUTe  ~  V2-poCT  ,  (11) 

where  is  the  peak  pressure,  is  air  density  (kg/m®),  c  is  the  speed  of 

sound  (m/s)  and  I  is  the  power  density  (W/m®).  After  recalling  the  relationship 
between  peak  and  average  pressure 


p  =  ^Peak  (-21 

and  inserting  this  into  Equation  (11)  and  solving  for  i  we  get 

p  ^ 

j  =  .  (13) 

PoC 


These  power  densities,  l,  were  then  integrated  over  the  surface  of  the  acoustic 
beam's  main  lobe  to  find  the  total  power,  P,^^,  in  the  lobe,  as  follows: 
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Expressing  this  integral  in  polar  coordinates 

-rp  dp  de 

Since  I  is  axially  symmetric,  the  Integral  becomes 

Ptot.i  =  2iifJ’"lpd9  .  (16) 

Simpson's  method  of  numerical  integration  evaluates  this  integral  as 

f^^flx)dx»S‘j(y^*4y^*2y2*4y^*  -  +2y„.2*4y„.i+y„)  .  (17) 

The  acoustic  power,  was  then  compared  to  the  electrical  Input 
power  applied  to  the  array,  The  transmitting  efficiency,  Ey,  was  the  ratio 
of  the  total  acoustic  input  power  divided  by  the  electrical  input  power  to  the 
array 

•  (1«) 

“input 


(14) 


(15) 


3.  Altitude  Dependence 

The  intensity  of  sound  in  the  atmosphere  depends  on  the  impedance 
of  air  (pc)  [Ref.  1 6],  where  p  is  the  air  density  (kg/m^)  and  c  is  the  speed  of 
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sound  (m/s).  These  variables  change  with  altitude,  since  p  depends  on  pressure 
and  temperature,  c  depends  only  on  temperature,  and  both  the  temperature  and 
pressure  vary  with  altitude.  At  higher  altitudes  and  lower  temperatures,  the 
array  transducers  are  less  efficient  than  they  are  at  sea  level.  This  does  not 
present  a  problem  if  the  echosounder  calibration  and  its  subsequent  use  in 
measurements  are  conducted  at  the  same  density.  However,  if  calibration  and 
operation  are  done  at  different  altitudes,  the  results  will  not  be  accurate.  The 
existing  acoustic  sounder  program  does  not  compensate  for  the  change  in 
transducer  efficiency  with  atmospheric  impedance,  pc. 

Figure  8  shows  the  experimental  equipment  arrangement  used  to 
determine  the  echosounder's  altitude  dependence.  A  signal  pulse  from  a 
HP3314A  Function  Generator  excited  a  Motorola  model  #  KSN  1005 A  ceramic 
piezoelectric  driver  which  was  sealed  inside  a  3*inch  diameter,  82  cm  long 
section  of  PVC  pipe.  The  same  3-inch  driver  received  the  pulsed  signal 
reflected  by  the  flat  surface  at  the  opposite  end  of  the  tube.  A  Nicolet  Pro  30 
Digital  Analyzer  measured  the  signal  intensity.  Measurements  were  made  as 
a  vacuum  pump  removed  the  air  from  the  pipe  over  a  pressure  range  of 
approximately  10-1000  mbar. 

A  plot  of  the  received  signal  voltage  versus  pressure  should  be  linear 
if  the  signal  amplitude  is  proportional  to  pc.  This  relationship  provides  an 
altitude  correction  to  the  anechoic  calibration.  The  AEOS  Hawaii  observatory 
at  3  km  altitude  has  a  pressure  that  is  30%  lower  than  at  sea  level. 
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4.  Modified  Echosounder  Equation 


Assume  that  an  echosounder  transducer  produces  a  peak  on-axis 
intensity,  at  a  distance  r,  given  by 

J„{i2)2(-^)  .  (18) 

Here  we  assume  that  the  far-field  intensity,  lo(ro)  was  measured  in  the 
laboratory  at  a  known  distance  r,,.  Using  the  far  field  amplitude  antenna  pattern 
F(9,0),  the  normalized  intensity  is  then 

|F{e,4>)  P  .  <«) 

where 

|F(0,0)P  =  1  (20) 

Consequently,  the  far-field  irradiance  on  a  target  at  range  r  will  be 

.  (21) 

•c  m 

If  the  atmospheric  volume  backscatter  coefficient  at  some  angle  f 
with  respect  to  the  incident  sound  direction  is 

o(C)  {m  steradian'^)  ,  (22) 

the  power  scattered  by  a  is 
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dP^(C)  =ItO(C)  dV  (W/sr)  ,  (23) 

where  the  ensonified  volume  element  of  area  da,  for  a  pulse  length  r  and  speed 
of  sound  c  is 

dv=  da  (m^)  .  (24) 

The  scattered  power  is  then 

=  Jo{^)2|F(e,<|))|2o(C)  (V/sr)  . 

X  M 

A  receiver  antenna  of  area  A  will  subtend  a  solid  angle  of  A/r’  from 
the  scattering  volume.  So  the  backscattered  power  received  at  f=/r/2  also 
depends  on  |F(  O,0)\^  and  will  be 

dp,(e,<i>)  ^da-4|F(e,4))i2  m 

(26) 

Using  dQ=da/r^  and  regrouping  the  antenna  pattern  terms,  the  backscattered 
power  becomes 

dp^{e,<l»  =  Jo  Vo{«/2) ^  |F(e,<l>)  l«dQ  (ko  .  (27) 

Dividing  by  A,  the  area  of  the  antenna,  we  can  rewrite  this  in  terms  of  the 
irradiance  at  the  receiver 
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(28) 


dj,(e,4>)  =Xo-^o(*/2) |F(e,<|))|"ciQ  {JL) 

Integrating  over  the  solid  angle  illuminated  by  the  antenna  and  including  the 
round  trip  absorption  loss  for  an  attenuation  coefficient,  a,  gives 
2 

J^  =  Jo  ^0(11/2)  ^e-2«f°^lF(e,4>)l"dQ  {-^)  .  (29) 

2  Jo  m2 

This  formulation  allows  us  to  use  a  laboratory  calibration  of  the  receiver 
sensitivity  expressed  in  terms  of  mV/Pa,  which  will  then  be  converted  to  Watts 
using  20x10®  Pa  -*  IxIO’^W.  The  solid  angle  covered  by  the  integral  should 
include  only  the  main  lobe,  since  the  array  enclosure  suppresses  the  sidelobes. 

Compare  this  with  the  previous  NOAA  formulation  of  Little,  Hall  or 

Neff 


^  =  [P^Er]  lo^  [R,f)]  [£1]  G]  .  (30) 

Although  similar  to  the  expression  above,  there  are  significant  differences.  Here 
the  terms  P„  E,  and  E,  as  well  as  G  are  virtually  impossible  to  determine  since 
they  depend  on  the  antenna  beamwidth  as  well  as  the  transducer  conversion 
efficiencies.  The  G  term  is  not  the  antenna  gain  but  is  related  to  the  F(P,0) 
integral.  This  equation  uses  the  total  power  P^  a  difficult  to  determine 
parameter  compared  to  the  peak  on  axis  irradiance. 
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IV.  EXPERIMENTAL  RESULTS 


A.  RECEIVING  EFFICIENCY 

A  B&K  Type  4228  pistonphone  calibrated  the  1/4”,  B&K  calibrated 
microphor^e  (serial  #1181 490) .  The  factory  calibration  of  the  pistonphone  was 
124.15  dB  for  a  reference  signal  of  20 //Pa.  Converting  acoustic  intensity  to 
pressure  gives 

^  124.1S<1B^ 

[10  20  ]  [20 xl0-«  Pal  =  32. 2±. 2 


for  the  signal  produced  by  the  pistonphone.  The  microphone  produced  a  signal 
amplitude  of  103.7  ±  .5  mV,»(,g  using  the  pistonphone  and  an  HP  3561 A  Signal 
Analyzer.  The  calibration  coefficient,  was  calculated  by  dividing  the 

microphone  measurement  by  the  pressure  generated  by  the  pistonphone 


„  _  103. 7  mV 

pietoa  -  32.25Pa 


3.22±.Q2 


mV 

Pa 


The  open  circuit  amplitude  of  the  signal  received  by  the  array  was  found 
to  be  504,7 ±  .5  for  a  1 .0  V  signal  input  to  a  1-inch  dime  tweeter  5.25 

m  from  the  array.  The  array  was  then  replaced  with  the  calibrated  microphone. 
The  microphone  signal  amplitude  was  24.8 ±.5 

Equation  (5)  states  that  the  array  receiving  efficiency,  E^,  should  specify 
the  efficiency  for  conversion  of  acoustic  power  back  to  electrical  power.  This 
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is  a  deceptive,  and  probably  erroneous  concept  since  the  array  is  a  transducer 
that  produces  a  voltage  from  an  applied  pressure.  The  21.1  ohm  impedance 
array  feeds  a  high  impedance  amplifier  <1  MQ).  The  amplifier  receives  very  little 
power  since  the  electrical  impedances  are  not  matched.  One  is  tempted  to 
compute  the  electrical  power  produced  from  the  acoustic  power  with 


r>  SMS 

~~Z~ 


where  Z  =  21.1  ohms  is  the  electrical  impedance  of  the  array,  but  this  is 
incorrect.  Instead,  the  receiving  efficiency  must  be  specified  in  terms  of  the 
open  circuit  voltage  produced  per  unit  pressure,  V/Pa,  as  is  done  for  calibrating 
microphones.  The  array  produced  a  signal  of  24.8 ±.5  Since  the 

microphone  calibration  factor  was  3.22  mV/Pa,  the  array  receiving  sensitivity 
was 


-(3.22^)  =65.5  ±1.3  (-^) 
24.8  Pa  Pa 


After  measuring  the  received  voltage  this  constant  allows  us  to  compute  the 
equivalent  pressure  and  then  the  power  using  the  relationship  20x10'°  Pa  -» 
1x10'’2  w/m^. 


B.  TRANSMITTING  EPFICIENCY 

The  1/4  inch,  B&K  calibrated  condenser  microphone  (serial  # 
1181490)  was  calibrated  with  the  pistonphone  as  before,  in  the  previous 
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section,  but  for  a  slightly  lower  ambient  atmospheric  pressure  (1013.0  mbar). 
The  pistonphone  factory  calibration  was  1 24.08  dB  for  a  reference  signal  of  20 
//Pa.  Converting  this  intensity  to  pressure  gives 

^  124.08  dB) 

[10  20  ]  [20xlO-^Pa^]  =31. 99 

the  signal  produced  by  the  pistonphone.  The  HP  3561 A  Signal  Analyzer 
measured  a  signal  amplitude  of  97.95  mV^Ms  when  the  microphone  was 
inserted  in  the  pistonphone.  The  calibration  coefficient,  was  3.062 

mV/Pa. 

Figure  9  shows  the  two  axes  of  symmetry  for  the  hexagonal  array.  The 
array  characteristics  were  measured  along  both  axes  of  symmetry.  TABLE  I 
shows  the  signal  amplitudes  measured  by  the  calibrated  microphone.  Four 
measurements  were  recorded  at  each  position,  two  side-mounted  and  two 
apex-mounted. 

The  measured  voltages  were  then  converted  into  power  densities  using 
Equations  (10)  through  (13).  The  total  power  contained  in  the  main  lobe  was 
calculated  using  Equations  (14)  through  (17).  TABLE  II  shows  the  measured 
power  in  the  main  lobe  for  each  array  orientation.  The  measurements  for  the 
side  mounted  array  were  greater  than  the  apex  mounted  because  of  the 
difficulty  in  aligning  the  array  maximum  peak  with  the  microphone.  The 
average  total  lobe  power,  P,^,^,  was  (3.0±  .3)x10'*  W. 
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Figure  9.  Power  density  measurements  of  the  hexagonal  array  were  made 
along  both  axes  of  symmetry. 


TABLE  I.  SIGNAL  AMPLITUDES  MEASURED  AT  DISTANCES 
INDICATED  FROM  MAIN  LOBE  CENTROID. 
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TABLE  II.  TOTAL  POWER  CALCULATED  IN  MAIN  LOBE  FOR  EACH 
ARRAY  ORIENTATION. 


1  Array  Orientation  | 

1  SIDE  Mounted 

1  Total  Power  (W) 

APEX  Mounted  | 

Total  Power  (W)  | 

Run  #1 

Run  #2 

Run  #3 

Run  #4  1 

0.000328 

0.000323 

0.000270 

0.000273  1 

The  input  power  applied  to  the  array  is  a  product  of  the  RMS  current  and 
voltage  applied  to  the  array  and  the  cosine  of  the  phase  angle.  An  HP4194A 
Complex  Impedance  Analyzer  found  the  phase  angle  to  be  0  =  -46.8  degrees 
at  5000  Hz.  Since  the  input  current  was  0.010  A„^  for  a  voltage  of  208.8 
mV^s,  the  array  input  power  Pi„p^  was  (1.43±.01)x10^  W.  The  transmitting 
efficiency,  Ej,  was  21  ±2  percent.  This  is  significantly  lower  than  the  50 
percent  that  Weingartner  [Ref  2]  estimated  and  was  used  in  the  echosounder 
program.  Substituting  this  value  for  E^  into  Equation  (71  would  increase  the 
recorded  C„^  values  by  a  factor  of  2. 

C.  ALTITUDE  DEPENDENCE 

The  HP3314A  function  generator  generated  a  1  volt,  5  kHz,  10  cycle 
signal  burst  which  was  both  transmitted  and  received  by  the  same  3-inch  driver 
inside  a  sealed  section  of  PVC  pipe.  TABLE  III  shows  the  sound  intensities 
received  by  the  3-inch  driver  as  the  air  was  pumped  out  of  the  sealed  PVC  pipe. 
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The  sound  intensities  were  plotted  versus  air  pressure.  Figure  10  shows  the 
sound  intensity  dependence  on  pressure  with  a  polynomial  best  fit 
superimposed  over  the  data  plot. 

Sound  pressure  readings  can  be  extracted  from  Figure  10  when  the 
applicable  air  pressure  is  known.  The  air  pressure  at  sea  level  is  approximately 
1013  mbar.  At  the  Starfire  Optical  Range  (SOR)  it  is  approximately  810  mbar, 
while  the  Advanced  Electro-Optic  Site  (AEOS)  has  an  air  pressure  of 
approximately  700  mbar.  These  pressures  translate  to  sound  voltages  received 
of  44.84  mV  (sea  level),  42.96  mV  (SOR)  and  40.74  mV  (AEOS)  for  an 
identically  transmitted  1 .0  V,  5.0  kHz  signal.  This  pressure  dependence  must 
be  factored  into  the  echosounder  calibration  before  it  can  be  applied  to  the 
AEOS  data. 
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TABLE  III.  SIGNAL  INTENSITY  MEASUREMENTS  RELATIVE  TO  AIR 
PRESSURE 


1  Pressure 

1  (mbar) 

Signal  Intensity 
(mV) 

Pressure 

(mbar) 

Signal  Intensity 
(mV) 

1  10.73 

471.28 

33.16 

1  14.12 

505.14 

34.8 

31.05 

0.9 

539.01 

36.06 

47.98 

2.04 

572.87 

36.8 

64.91 

3.18 

606.73 

37.8 

81.84 

4.62 

640.6 

39.1 

98.78 

6.06 

674.46 

40.01 

132.64 

9.46 

708.33 

41.06 

166.5 

12.6 

742.19 

41.65 

200.37 

16.0 

776.05 

42.52 

234.23 

19.06 

809.92 

42.94 

268.09 

21.5 

843.78 

43.3 

301 .96 

23.34 

877.64 

43.52 

1  335.82 

26.1 

911.51 

43.92  1 

1  369.69 

27.6 

945.37 

44.25 

1  403.55 

30.18 

979.24 

44.94 

1  437.41 

31.6 

1013.1 

44.84 

1  454.35 

32.38 
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Voltage  (mV) 


Figure  10.  Signal  intensity  dependence  on  air  pressure  with  polynomial  fit 
superimposed  on  data  plot. 
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D.  MODIFIED  ECHOSOUNDER  EQUATION 

Measurements  of  the  peak  sensitivity  of  the  array  were  very  sensitive  to 
alignment  of  the  array  since  small  angular  misalignments  moved  the  maximum 
away  from  the  measurement  point.  The  highest  values  in  TABLE  I  reveal  this 
phenomena.  Consequently,  the  higher  values  of  TABLE  I,  the  first  two 
columns,  are  probably  more  reliable  than  the  second  two  columns  that  have 
some  misalignment  error.  The  first  column  from  TABLE  I  shows  a  peak  on  axis 
microphone  voltage  of  1.37  mV  for  0.208  mV  applied  to  the  array.  This 
corresponds  to  6.49  mV  for  1 .0  volt  applied  to  the  array.  Dividing  by  the 
calibrated  microphone  sensitivity  of  3.062  mV/Pa,  this  corresponds  to  2.1 2  Pa 
or  1.12x10'^  W/m^  at  5.25  m  from  the  source.  Referencing  this  to  1  m  the 
transmitting  efficiency  becomes 

=  0-309  { — 

Calculating  the  average  transmitting/receiving  antenna  angular  distribution 
over  the  main  beam  lobe  was  a  bit  more  involved.  The  data  of  TABLE  I  were 
normalized  so  that  F(0,0)  =  1  and  dO  =  sin8  69  d0.  The  integral  was  simplified 
by  assuming  symmetry  in  the  <p  direction  and  using  By  the  angle  of  the  first 
diffraction  pattern  minimum, 
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|F(e,4>)  I*  dO  =  I F(e)  lysine ded(>  =  211  J*^®'lF(e)  l^sinOdD  . 


Using  Simpson's  method  the  integral  was  evaluated  numerically  and  found  to 
be  0.0096 ±.0004. 

Recalling  the  transmitting  efficiency  of  0.309  W/|Vm)^  (at  1  m  and  1  Vr^s) 
and  the  angular  integral  of  0.0096  (rounded  off  to  0.01 )  allows  us  to  compare 
the  previously  stated  Equation  (29) 

■rr  =  -ro-^0(*/2)  -^e-2«^f“^|F(e,<fr)|"cfQ  (-i^)  ,  (29) 

2  Jo  np- 

with  the  older  calibration  expression.  Equation  (30) 


=  [e-2«*]  [o^(F,f)}  [^]  I-4g1 

2  J?2 


(30) 


Dividing  Equation  (30)  by  the  antenna  area.  A,  and  eliminating  like  terms  in 
both  equations,  allows  us  to  compare  the  sound  Intensities,  I,,  at  r  =  1m 


I±. 

AE^ 


—  PjEj, G  . 


Substituting  values  for  V^s»  (25  V),  Et  (0.5)  and  G  (0.4)  we  have 


(25  V)^ 
21  0 


•  (0.5)  •  (0.4)  =6.0  (J^) 

nr 
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for  the  transmitted  sound  intensity.  This  intensity  was  multiplied  by  the 
receiving  efficiency,  E^,  array  area,  phase  angle  and  the  impedance  (21  Q)  to 
find  the  received 

6.0  (-^)  -  (0.5)  •  (0.0866  7172)  -cos  {-46. 8®)  •  (21  O)  =  3 . 76  ( 
nr 

for  comparision  with  the  modified  echosounder  equation. 

For  Equation  (29),  using  r„  =  1  m,  we  have 

=  Ctrsne  '  W  <())  |‘  dQ 

Substituting  values  for  (0.309),  (25  V)  and  the  angular  integral 

(0.01 ),  we  have 

0.309  (<*y  (W77)2)  '  (25  V)2  -0.01  =  1.93  (  —  ) 


for  the  transmitting  sound  intensity.  A  receiving  conversion  coefficient  was 
found  using  the  receiving  sensitivity  (65.5  mV/Pa),  mentioned  earlier,  and  the 
pressure-to-intensity  conversion  i.e.. 


,  6.55x10-2  V.  2  _  /  y2 

20X10-*  '  ■  W/m^ 


) 


The  transmitting  voltage  was  multiplied  by  the  receiving  conversion  coefficient, 
the  area  of  the  array  and  the  phase  angle  to  give  the  receiving  for 
comparison 
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1.93  {-^)  -1.72  (—^)  =3.32  (1^2) 

W/m^ 

Comparing  these  two  expressions  we  find  that  they  are  approximately  equal. 
Since  both  of  these  approaches  result  in  the  same  approximate  value,  it  is  not 
clear  why  a  factor  of  four  discrepancy  between  the  echosounder  and 
microthermal  probe  values  exists. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

To  date,  most  acoustic  echosounders  have  used  an  expression  for 
computing  Ct^  and  optical  turbulence  that  contain  transducer  total  power  and 
efficiency  terms  Pj,  E^,  E^,  and  G  that  depend  on  the  transducer  electrical- 
acoustic  efficiencies  as  well  as  the  antenna  beam  pattern  coupled  together. 
This  coupling  has  made  calibration  exceedingly  difficult.  To  avoid  this  problem, 
this  thesis  research  considered  two  separate  calibration  approaches.  The  first 
was  to  measure  C/  directly  using  microthermal  probes  on  a  tower  and  compare 
these  results  with  the  acoustic  sounder  measurements.  The  second  approach 
was  to  recast  the  calibration  problem  into  separate  transducer  transmitting  and 
receiving  calibration  coefficients  measured  in  an  anechoic  chamber.  These 
calibrations  include  the  work  of  Probert-Jones  that  state  that  a  monostatic 
radar  must  include  an  average  over  the  product  of  the  transmitting  and 
receiving  antenna  beam  patterns  leading  to  an  integral  over  the  forth  power  of 
the  antenna  amplitude  angular  distribution. 

The  September  1992  tower  measurements  showed  that  the  original 
acoustic  sounder  calibration  expression  overestimated  the  actual  atmospheric 
turbulence  by  a  factor  of  four.  This  factor  of  four  had  been  seen  in  previous 
work,  but  the  source  of  this  factor  was  not  clear. 
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Using  the  acoustic  calibration  procedure  in  an  anechoic  chamber  the 
transmitting  coefficient  was  0.309  W/(Vm^)  at  1  m  for  each  volt  RMS  applied 
to  the  array.  The  receiving  sensitivity  was  65.5  ±  1 .3  mV/Pa.  Combining  these 
two  coefficients  produced  results  that  were  nearly  identical  with  the  previous 
calibration  expression.  To  date  these  coefficients  have  failed  to  resolve  the 
factor  of  four  discrepancy  between  the  acoustic  calibration  measurements  and 
calibration  measurements  performed  using  microthermal  probes. 

This  thesis  also  investigated  the  dependence  of  the  acoustic  echosounder 
calibration  on  the  acoustic  impedance  of  the  air,  pc,  which  varies  with 
atmospheric  temperature  and  pressure.  Although  a  pc  decrease  of  30  percent 
was  expected  from  sea  level  to  3  km,  the  actual  driver  measurements  found  a 
9  percent  reduction.  These  occurred  since  the  acoustic  drivers  operated  at 
resonance  and  the  displacement  of  air  around  the  driver  horn  phasing  plug 
provided  considerable  damping.  This  damping  declined  with  reduced  air 
pressure  leaving  the  driver  displacement  relatively  unchanged. 

We  consider  the  tower  calibrations  to  be  the  most  accurate  at  this  time 
since  they  are  consistent  with  previous  turbulence  measurements.  The 
inconsistency  between  the  microthermal  probe  and  acoustic  calibrations 
remains  intriguing. 
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B.  RECOMMENDATIONS 

The  factor  of  four  discrepancy  between  the  acoustic  and  microthermal 
echosounder  calibrations  needs  to  be  resolved.  The  source  of  this  could  be  as 
simple  as  (an  elusive)  programming  error.  Although  this  has  been  examined, 
carefully,  the  potential  for  such  an  error  exists.  On  the  other  hand,  this 
discrepancy  may  indicate  a  more  fundamental  problem  with  the  theoretical 
expressions  relating  atmospheric  turbulence  to  the  acoustic  scattering  cross- 
section  or  that  atmospheric  turbulence  in  the  3  cm  size  range  is  four  times 
larger  than  expected.  Another  possibility  is  that  the  micothermal  reference 
probes  are  in  error  by  a  factor  of  four. 
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